Evolution and Symmetry of Multipartite Entanglement 
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We discover a simple factorization law describing how multipartite entanglement of a composite 
quantum system evolves when one of the subsystems undergoes an arbitrary physical process. This 
multipartite entanglement decay is determined uniquely by a single factor we call the entanglement 
resilience factor (ERF). Since the ERF is a function of the quantum channel alone, we find that 
multipartite entanglement evolves in exactly the same way as bipartite (two qudits) entanglement. 
For the two qubits case, our factorization law reduces to the main result of Nature Physics 4, 99 
(2008). In addition, for a permutation P, we provide an operational definition of P-asymmetry of 
entanglement, and find the conditions when a permuted version of a state can be achieved by local 
means. 
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With the emergence of quantum information science 
in recent years, much effort has been given to the study 
of entanglement [H, 0- It was realized that highly en- 
tangled states are the most desirable resources for a va- 
riety of quantum information processing (QIP) tasks, 
such as quantum teleportation Q , super-dense coding [J] , 
entanglement-based quantum cryptography [5| , error cor- 
recting codes J6J , and more recently, one-way quantum 
computation [7|. Due to the effect of decoherance in- 
duced by the coupling of the subsystems with the en- 
vironment, the entanglement of the composite quantum 
system decreases in time. It is therefore critical, for the 
implementations of many important QIP tasks, to under- 
stand the behaviour of entanglement under the influence 
of decoherance or noise. 

To study the evolution of entanglement it seems to 
be necessary first to study the evolution of the quantum 
state describing the composite system and then to cal- 
culate its entanglement. For example, a situation where 
no energy is exchanged with the environment, the mas- 
ter equation involving the Lindbland operators can be 
used to determine the state evolution. Indeed, the elab- 
orate theory on state evolution was the method used 
by many researchers (e.g. see references in Q). How- 
ever, the drawback of this technique is that for mul- 
tipartite systems (or higher dimensional systems) the 
state equation can be very hard to solve and therefore 
the evolution of entanglement can be determined only 
in very special cases. Quite recently, a new way was 
found 0] (see also [j| for an earlier similar work) to char- 
acterize the evolution of entanglement in two qubits sys- 
tems, by which the evolution of concurrence [10( (a two 
qubit measure of entanglement) is determined directly in 
terms of the evolution of a maximally entangled state; 
i.e. a Bell state. This technique was generalized to de- 
termine the evolution of the G-concurrence [ll| of two 
qudits in [lii |. In both Q and [IH the authors used the 
Choi-Jamiolkowski isomorphism in order to derive the 



equations describing the time evolution of entanglement. 
Hence, since the Choi-Jamiolkowski isomorphism applies 
only for bipartite systems, it may give the impression 
that such entanglement-evolution equations can not be 
extended to multipartite settings. 

In this Letter we discover a simple factorization law 
describing how multipartite entanglement of a compos- 
ite quantum system evolves when one of the subsystems 
undergoes an arbitrary physical process. Quite remark- 
able, this factorization law holds for arbitrary number 
of parties, and reduces to the factorization law given 



12] for the bipartite case. Our key idea is to use 



measures of entanglement that are invariant under the 
group G = SL(di,C)®SL(d2,C)®---®SL(d n ,C), where 
d\,d,2, d n are the dimensions of the n-subsytems, and 
SL(d, C) is the group of d x d complex matrices with 
determinant 1. The group G represents (determinant 
1) stochastic operations assisted by classical commu- 
nications (SLOCC) and has been used extensively in 
the classifications of multi-partite entanglement. It is 
therefore clear from our analysis that even in the bipar- 
tite case, it is the invariance under G, rather than the 
Cho-Jemiolkowski isomorphism, that is necessary for the 
derivation of the factorization law. 

In addition to the factorization law, we also provide an 
operational definition of P-asymmetry of entanglement: 
a multipartite entangled state contains P-asymmetric en- 
tanglement if its subsystems can not be permuted (ac- 
cording to the permutation P) by means of LOCC. We 
show that in general states have P-asymmetric entangle- 
ment, and by using measures of entanglement that are 
invariant under G, we are able to generalize the main 
result of [HI to the case of multi-partite systems. 

While two-party entanglement was very well studied, 
entanglement in multi-party systems is far less under- 
stood. Perhaps one of the reasons is that n qubits (with 
n > 3) can be entangled in an uncountable number of 
ways PJ, LLa ll2| with respect to SLOCC. It is there- 
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fore not very clear what role entanglement monotones 
can play in multi-qubits system unless they are defined 
operationally. One exception from this conclusion are en- 
tanglement monotones that are defined in terms of SL- 
invariant polynomials IB- 23 1. 

For example, for two qubits the concurrence [l(| is 
an entanglement monotone that is invariant under the 
action of the group SL(2,C) O SL(2,C). That is, for 
a given two qubits state \ip) the concurrence of \ip) is 
the same as the concurrence of the unnormalized state 
A ® B\ip), where A and B are 2x2 complex matrices 
with determinant 1. The concurrence originally was de- 
fined as a mathematical tool to calculate an operational 
measure of entanglement, namely, the entanglement of 
formation. Nevertheless, the concurrence, as being the 
only SL(2, C) <g> SL(2, C) invariant 2-qubits homogeneous 
measure of degree 1 , was used as the key measure of en- 
tanglement in numerous QIP tasks (e.g. see [l], Q and 
references therein). 

Another example is the square root of the 3-tangle 
(SRT) [H. The SRT is the only SL(2,C) <E> SL(2,C) (g) 
SL(2, C) invariant measure of entanglement that is ho- 
mogenous of degree 1. The SRT also plays an important 
role in relation to monogamy of entanglement, and in 
fact it is the only measure that capture the 3-way entan- 
glement. For 4-qubits or more, the picture is different 
since there are many homogenous SL-invariant measures 
of entanglement, such as the the square root 4-tangle [3] 
or the 24th root of the Hyperdeterminant [22|. We now 
define all such measures that will be satisfying our fac- 
torization law; such measures were first discussed in 17 1 . 



Definition 1. Set G = SL(di,C) <g> SL(d 2 ,C) ® ••• <g> 
SL(d„,C), U n = <C dl (g>C d2 ® - • -®C<S and B{H n ) the set 
of all bounded operators (e.g. density matrices) acting 
on % n . A SL-invariant multi-partite measure of entan- 
glement, Ei nv , is a non-zero function from B(U n ) to the 
non-negative real numbers satisfying the following: 

(i) It is G-invariant; that is, 

Ei n v(gpg*) = E inv (p) , 

for all g e G and p G B^H n ). 

(ii) It is homogenous of degree 1; i.e. 

Einv {rp) = rE inv (p) 

for all non- negative r and all p e B(H n )- 

(iii) On mixed states it is given in terms of the convex 
roof extension; That is, 



E n 



,(ipi 



where the minimum is taken with respect to all pure 
states decompositions of p = "YlnPiipi (here ipi = {tpi\ 
is a rank 1 density matrix). 



Remark. The criteria in the definition above guarantee 
that Ei nv is an entanglement monotone 17J. Note also 
that the construction via the convex roof extension is con- 
sistent with conditions (i) and (ii). The concurrence, the 
G-concurrence, and the SRT are all satisfying the con- 
ditions in the definition above. It can be easily checked 
that Ei nv is unique (up to multiplication by a positive 
constant) for the bipartite case with d\ = c?2 and for 
three qubits, but it is not unique for n-qubits with n > 3. 
Indeed, for 4 qubits there are 4 algebraically independent 
SL-invariant polynomials that generate a whole family of 
such SL-invariant measures [la, [l9l - l2ll ]. Note however 
that if the dimensions of the subsystems {di} are not all 
equal then a SL-invariant measure E inv may not exists. 
For example, in the bipartite case with di ^ d 2 , E inv 
does not exists. For three parties, on the other hand, 
with d\ = di = 2 and dz — 3, such a measure exists; it is 
given in terms of the Hyperdeterminant 22j . 

The criteria in the definition above are motivated from 
two observations. First, if two states are connected by 
SLOCC reversible transformation, then essentially they 
both must have the same type (though not necessarily 
the same amount) of multi-qubit entanglement. That is, 
let \tp) £ U n and let \tp) = g\ip)/\\g\i())\\, where g € G. 
Then, criteria (i) and (ii) guarantee that if Ei nv (\ip)) ^ 
then also Ei nv (\<p)) ^ 0. The second observation says 
something about the amounts of the multipartite entan- 
glement in l^) and \<p). It is based on some results first 
discovered in [13] and discussed further in [r|. We now 
describe shortly these results. 

Let \ijj) £ U-n and consider the set of (in general 
non-normalized) states G\ip) (i.e. the orbit of \ip) un- 
der G). By definition, if is generic then G\ip) is 
closed. Therefore, for most states G\if>) is closed. If 
G\tj}) is not closed then consider its closure, and de- 
note by \(p) the state in G\tp) with the minimum norm; 
that is (<p\f> ) < ( <p\<p) for all \4>) <G G\ip). The state 
\ip) = \<p)/y/ (<p\tp) is called a normal form (see also 
the critical set in appendix A of [3]). Moreover, note 
that if the normalized state \tp) is a normal form, then 
HfflV'MI — 1 f° r au 9 ^ G, with equality if and only if 
g e SU(di ) ® SU(d 2 ) ® ■ • • ® SU(d„) 0, E3- 

The properties (i) and (ii) in the definition above imply 
that if \tp) is a normal form, and if \ip) — g\ip)/\\g\ip)\\, 
then 



E inv (\(p)) 



Einv(g\lp}) E mv (\ip)) 



\g\m 



< E lnv {\ip)) 



That is, criteria (i) and (ii) imply that among all the 
states that can be obtained from \tp) by SLOCC, the nor- 
mal form \ip) has the maximum amount of Ei nv - Indeed, 
in [l6|; [ItJ it has been shown that \ip) is a normal form 
if and only if each qudit is maximally entangled with the 
rest of the qudits (i.e. the local density matrices of all 
qudits are proportional to the identity). Therefore, crite- 
ria (i) and (ii) are consistent with this result, and we can 
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consider the normal forms as maximally entangled states. 
More details and further motivation for the first criterion 
can be found in the extensive literature on the character- 
ization of entanglement in terms of SL-invariant polyno- 
mials (see for example flil - l22| and references therein). 
Our last remark on Def. [T]we summarize in the following 
lemma. 

Lemma 1. Let \ip) 6 ~H n and E inv as defined in Def. [7J 
Then, for a matrix M : H n — > H n of the form M — 
A 1 ® A 2 ® ■ ■ ■ ® A n we have E in0 (M\ij)) {ip\M^) = if 
there exists 1 < k < n such that det = 0. 

Proof. Without loss of generality, assume det^i = 0. 
Denote by |0) a normalized vector in C rfl such that 
Ai|0) = 0. Denote by \k) (with k = 1, 2, d ± - 1) other 
vectors in C dl , completing |0) to a basis. With this basis 
we can write 

M\>l>) = £ \ v k)\<Pk) , with I^gC*®-.-®^, (1) 

k=l 

and \vk) = Ai\k). Note that \vk) and \(fk) are not nec- 
essarily normalized and that the sum above starts from 
k = 1 and not from k = since Ai|0) = 0. Denote by P r 
the projection to the space V = sp&n{\vk)} dl =1 1 , where 
r = dim V. For < t £ R, let D t be the following di x di 
matrix: 

D t = t r/{r - dl) (I - P r ) + tP r . 

Denote also g t = D t <g> I <g> ■ • • <g> I. Clearly det g t = 1; i.e. 
g t £ G. Hence, 

E inv [M\i>)(iP\M*] =E inv \g t M\i))^\M^g\ 
= E inv [t 2 M\i)){^\M^] = t 2 E inv [M|V)(^|Mt] , 

where we have used Eq.([T]) and criteria (i) and (ii) of 
Def. Q] Since the above equality is true for all t > 0, we 
must have E mv [M\tp) (ip\M^] =0. □ 

We are now ready to discuss the main result of this 
paper. 



Definition 2. Let $ : B(C d ) -)• B(C d ), be a quantum 
channel acting on d x d positive semi-definite matrices 
(i.e. density matrices). Any such channel has Kraus 
representation $(■) = J^j with Kraus operators 

^2jKjKj < 1. We define the entanglement resilience 
factor of $ to be 



Note that < .F[$] < 1 due to the geometric- 
arithmetic inequality and the fact that for all Kraus rep- 
resentations of $, ^2jKjKj < 1. Recall also that all 
Kraus representations of a quantum channel are related 
by a unitary matrix. In the theorem below we give an 
operational interpretation for J-"[$|. 

Theorem 2. Let \ip) £ H n and \<j>) £ H n be two states 
with non zero value of Ei nv . Denote by A = $ <£) I <g> • • • ® 
I, where $ is an arbitrary quantum channel, which may 
represent the influence of the environment on the first 
audit. Then, 



E mv [A (|V>)M)] 

Einv(\ll>)(ll>\) 



E inv [A(\<l>) (</>{)] 



= J-[$]. (3) 



That is, the ratio between the final entanglement and the 
initial entanglement depends solely on the entanglement 
resilience factor of the channel. 

Remark. In the bipartite case with d\ = d 2 the formula 
above reduce to the one given in [g, 12 1, by replacing the 
state with a maximally entangled state, and by taking 
Ei m) to be the concurrence II Oil or G-concurrence 




two qubits or two qudits, respectively. However, in 
the entanglement resilience factor (ERF) of the channel 
was not introduced. A remarkable observation is that 
the ERF depends only on a single qudit channel. Now, 
consider the case of two qubits. Then, by taking Ei nv to 
be the concurrence C and replacing above with a Bell 
state we get the following formula for the ERF: 



^[$] = c($®i(|^ + )(^ + |)) 



(4) 



which can be determined completely by using the Woot- 
ters formula. It is remarkable that for any number of 
qubits and for any choice of E inv , this is the unique for- 
mula that is needed to be calculated in order to determine 
the ERF of a channel acting on a qubit. Similarly, for a 
channel $ acting on a qudit, the ERF is given in terms 
of the G-concurrence [llj : 



^[$] = G($®I(|^+)(V' + |)) 



(5) 



where here stands for a maximally entangled state 
in C d <g> C d . Eq.© provide the unique value of the ERF 
and can be used to the determine the ratios in Eq.Q 
independent on the choice of Ei nv or the number of qudits 
involved. 

Proof. Given the channel $(■) = Kj(-)K^ , the den- 



sity matrix p = ($ <g) I <g) • • 
pure state decomposition: 



<8 I) |V')(V'I have the following 



= min ^ | det if. 



,2/d 



(2) 



P = J2\vj){vj\ , |%) sit,-® I® •••»%) 



where the minimum is taken with respect to all the Kraus 
representations of $. 



Note that \vj) are not normalized. Moreover, denote by 
p = \wi)(wi\ the optimal decompositions of p. That 
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is, denote pi = (wi\wi) (so that p i \wi) is normalized) 



E inv (p) = ^ pjE, 



mr I —\Wi)(Wi 

Pi 



Einv (\Wi)(Wi\) 



where we have used the fact that E inv is homogeneous of 
degree 1. Now, since and are two deferent 

decompositions of p, they are related to each other via 
a unitary matrix U. That is, if the two sets and 
{|ii)t)} do not have the same number of vectors we add 
zero vectors to the smaller set and then we have 

\wi) = Uij\vj) = Mi <g> I® • • • <g> I|^>) , 

3 

where U is a unitary matrix, and Mi = ■ UijKj form 
another Kraus representation to the same quantum chan- 
nel $. Now, w.l.o.g. (see lemma [TJ we can assume that 
detMj ^ 0. Hence, we can write 



\vh) = (det Mi) 1/rf 



M, 



(detMi) 1 /^ 



® I® ••• ®I 



Since Ei nv is G-invariant and homogeneous, we get 
E inv {\wi)) = \det Mi\ 2/d E inv (\ij))) and thus 



Einv(p) = E \tetMi\ 2/d E ino (\il>)) . 

i 

What is left to show is that 

=J2\ detM i\ 2/d ■ 



(6) 



(7) 



To see that, note that the unitary U has been chosen 
such that the decomposition p = J2i \wi)(wi\ is opti- 
mal. Based on Eq.([6j, U has been chosen such that 
Zi\detMi\ 2/d gets the minimum possible value among 
all the different Kraus representations of $. Hence, the 
equality in Eq.(0 must hold. □ 

The measure Ei nv is a convex function since it is de- 
fined in terms of the convex roof extension. Thus, we 
have the following corollary. 

Corollary 3. Let p G H n be a multipartite mixed state 
with non-zero value of Ei nv , and let A and $ be as in the 
theorem above. Then, 



E inv [A (p)} 

Einv (p) 



<-F[S] 



The following corollary is an immediate consequence 
of the equation above. 

Corollary 4. Let p G Ti. n be a multipartite mixed state 
with non-zero value of Ei nv , and let {$fc}fe=i,2,...,n be a 
set of n quantum channels. Then, 



Ej nv [$1 ® $ 2 8) ■ ■ ■ $n (p)} 
Einv (P) 



< 



(8) 



As a simple illustration of the above theorem and corol- 
laries, consider the case of three qubits. In three qubits 
the only G-invariant measure of entanglement, Ei nv , is 
given by the SRT [23| on pure states, and on mixed states 
it is defined in terms of the convex roof extension. The 
GHZ state maximize this measure. Applying the theorem 
above to this measure gives 



Ei; 



>I\GHZ)(GHZ\) = J"[$] , 



where J-[$] can be calculated via Wootters formula (see 
Eq.Q). That is, we have found a closed formula for 
the SRT for all mixed states of the form p = $ eg) I <£> 
I\GHZ)(GHZ\. Moreover, the corollaries above provides 
upper bounds for Ei nv on states of the form p = $i Cg) $2 ® 
% 3 \GHZ){GHZ\. 

One can also ask how the multipartite entanglement 
evolve after a separable measurement is performed by 
the n parties. In the following lemma we obtain an up- 
per bound on the ratio between the initial entanglement 
and the final average entanglement after such a measure- 
ment is performed. This lemma will be very useful in our 
discussion on symmetry of multipartite entanglement. 

Lemma 5. Let A(-) = 53fe-^fc(')-^fe be a trace preserving 



separable operation, where 
andY,k M t M * ^ L Then ' 



A 



(fc) 



(fc) 



^kPkEinvi&k) 
Einv (p) 



< 



E 



detMJ 2 / d 



(9) 



where = —M^pM, and pk = TrMkpMl. Further, 



k 



detM h U 



E 

k 



detA 



(fe)i 



detA«|3 < 1 



(10) 

with equality if and only if all the operators {A^} are 
proportional to unitaries. That is, if J\ | det Mk\ 2 / d = 1 
then A is a mixture of product unitary operations. 

Proof. The upper bound in Eq. {[9j) is a direct conse- 
quence of conditions (i) and (ii) in Def. [TJ and the upper 
bound in Eq. (|10[) is a direct consequence of the geometric- 
arithmetic inequality. □ 

We now show that SL-invariant measures of multipar- 
tite entanglement can also be very useful to determine 
the symmetry of multipartite entanglement. 

Definition 3. Let P be a permutation on n parties. 
Let p : Tin — > "Hn be a multipartite density matrix and 
denote by VppVp the "permuted version" of p, where Vp 
is unitary operator that permutes the subsystems. Then, 
the entanglement contained in a multipartite state, p, is 
said to be P-Symmetric, if by LOCC we can produce the 
permuted version of the state; i.e. VppVp. 



k=l 



The theorem below generalize the main result of 
to multipartite states. 



13] 
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Theorem 6. Let Ei nv in Def. \1\ be also invariant under 
some permutation P of the n qudits. Let p : % n — > % n be 
a mixed state for which E inv (p) > and assume that the 
entanglement of p is P-Symmetric. Then, the permuted 
version of the state can be achieved by some product uni- 
tary operation Uai ® Ua 2 ® ■ • • <8> Ua„ ■ 

Note in particular that if a state with E inv > has 
different entropies of subsystems, it can not be permuted 
by LOCC, as local unitaries can not change local en- 
tropy. However, if a state has distinct local entropies, 
but Ei nv = 0, then it may be possible to permute the 
state by LOCC. As a simplest example consider the state 
p = PA! <S> PA 2 ® • • • ® PA n ■ Clearly, any such state, or its 
permuted version, can be generated locally even though 
its local entropies can be distinct. Note also that in three 
qubits, the SRT is also invariant under permutations [23j . 
and therefore can be used for the theorem above for 
all permutations P. Same is true for all E inv that are 



also permutation invariant, such as the 4-tangle [18[ and 
all Hyperdeterminants 22 1 (see also [2(| for other such 
Einv)- The proof below of theorem [5] is based on lem ma p 
and follows the exact same lines as in Theorem 1 of (13| . 

Proof. Suppose E(p) > and let A (as defined in 
lemma [5]) represents an LOCC map that permutes p; i.e. 
A(p) = VppVp. Hence, J2kPk°~k = VppVp. Using the 
invariant of Ei nv under the permutation P we get 

Einv(p) = E inv {V P pVl) = E inv (^2p k a k ) 

k 

< ^PkE mv {a k ) < I det M k \ 2/d E mv (p) 

k k 

Since we assumed that Ei nv (p) > we get that 
Y, k I detM k \ 2/d > 1. Thus, from lemma [S] we conclude 
that A must be a mixture of product unitaries: 

Hp) = ^2pkU k pUl = ^2p k (Jk , 



where Uk = Ua} ® UaI ® ' ' ' ® ^aI IS a product of 
unitaries. Hence, all the states o~k have the same von 
Neumann entropy S as p, and we get 



•5 5>**k = S (y P pv£) = S(p) = Y,PkS(<J k ) • 

V k I k 

From the strict concavity of the von Neumann entropy 
we get that all the o~k are the same and therefore Vp — 



U 



(i) 



(1) 



• ® Ujp is a product of unitaries. 



□ 



In conclusions, Eqs. (|3I8[) provides us, for the first 
time, with closed expressions for the time evolution of 
multipartite entanglement of a composite system inter- 
acting locally with the environment. These expressions 
emerge from the SL-invariance of the measures defined 



in Def. Q] and not from the Jamiolkowski isomorphism 
which is the methodology used in Ref. 0, 12 1. Amaz- 
ingly, the evolution of multipartite entanglement (under 
one local channel) is determined completely by the ERF 
defined in Eq.@, irrespective to the number of qudits 
in the system. For multi-qubits systems, the ERF has 
a closed formula given in terms of Wootters concurrence 
formula. In other words, there is no need to solve any 
master equations in order to determine the time evolu- 
tion of multipartite entanglement. 
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